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ABSTRACT 

We present radio observations of the unique, recently formed, planetary nebula (PN) associ- 
ated with a very long -period OH/IR variable star V1018 Sco that is unequivocally still in its 
asymptotic giant branch phase. Two regions within the optical nebula are clearly detected in 
nonthermal radio continuum emission, with radio spectral indices comparable to those seen 
in colliding-wind Wolf-Rayet binaries. We suggest that these represent shocked interactions 
between the hot, fast stellar wind and the cold nebular shell that represents the PN's slow 
wind moving away from the central star. This same interface produces both synchrotron radio 
continuum and the optical PN emission. The fast wind is neither spherical in geometry nor 
aligned with any obvious optical or radio axis. We also report the detection of transient H2O 
maser emission in this nebula. 

Key words: masers - stars: AGB and post-AGB - planetary nebulae: general - radiation 
mechanisms: nonthermal 



1 INTRODUCTION 

Cohen, Parker & Chapman (2005: hereafter CPC) have reported 
the discovery of a very faint planetary nebula (PN) centred on the 
long-period pulsating OH/IR star OH 354.88-0.54, also known as 



V1018 Sco. The PN was found during examination of survey field 
HA630 of the Southern Ha Survey (SHS: Parker et al. 2005a) by 
QAP, as part of a systematic search for new Galactic Plane PNe 
(e.g. Parker et al. 2005b, 2006). Its OH characteristics mark this 
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star as an intermediate mass asymptotic giant branch (AGB) star, 
with initial mass > 4Mq . After rejecting the possibility that the 
PN arises from shock interaction with the surrounding external 
medium, or that the star is a "born-again" (e.g., Iben 1984) ob- 
ject that has undergone a very late thermal pulse (e.g., Blocker & 
Schonberner 1997) as it began to evolve along its cooling track, 
CPC speculated that this object represents a phase of PN evolution 
that has been formerly unobserved. They concluded that this star's 
fast wind has very recently turned on and it is the interaction be- 
tween this wind and the previously shed outer layers of the star that 
produces the ionization that reveals the optical PN. 

To explore this intriguing phenomenon further we have un- 
dertaken a series of radio continuum and water maser observations 
with the Australia Telescope Compact Array (ATCA), originally 
with the objective of mapping out the spatial distribution of ionized 
material within the PN. Pre-existing radio data are discussed in §2. 
§3 describes the new ATCA continuum and H2O maser observa- 
tions, while §4 shows an enhanced optical image of the PN. The 
discussion follows in §5, with our conclusions in §6. 

2 PREVIOUS RADIO CONTINUUM OBERVATIONS 

The best estimate for the position of the star itself is: il h 35 m 
02.73 s ±0.02 s , -33° 33' 29.41"±0.24" (CPC). Neither the PMN 
or NVSS surveys includes any continuum source within several ar- 
cmin of OH 354.88-0.5, while the VLA survey of the inner plane at 
5 GHz by Becker et al. (1994) is limited to |b|<0.4°. The 1.4-GHz 
continuum survey of small diameter sources by Zoonematkermani 
et al. (1990) does include a bright, compact radio object (2.4" in 
size) within 7" of the maser position. However, this survey com- 
prises two 3-MHz wide bands centred at 1441.5 and 1611.7 MHz, 
so the higher frequency band is contaminated by the strong 1612- 
MHz maser emission. This problem was recognized by White, 
Becker & Helfand (2005),who give an updated position at 20 cm 
of ll h 35 m 02.590 s -33° 33' 31.18" (J2000) and warn of this 
contamination of their band-averaged flux densities. Zijlstra et al. 
(1989) quote an upper limit of 1.4 mJy at 2 cm from a VLA obser- 
vation in September 1986. 

Consequently, the only radio survey to have correctly de- 



tected and assessed the continuum from OH 354.88-0.5 is the first 
Molonglo Galactic Plane Survey (MGPS1: Green et al. 1999), at 
843 MHz. Extracting the individual MGPS1 images, from which 
the final survey images were constructed, shows a source at the 
maser position. Fitting a 2D-Gaussian to the object indicates an 
unresolved source with an 843-MHz flux density of 10.8±0.6 mJy, 
located at ll h 35 m 01.52 s -33° 33' 36.90" (J2000) The beam size 
of the Molongolo Observatory Synthesis Telescope at this location 
is 75" (N-S) by 43" (E-W) so the radio emission could occupy the 
full extent of the Ha ring (roughly a 39" circle as shown by the 
images in §4) yet still be unresolved in MGPS1 images. Therefore, 
observations with the ATCA were planned. 

3 NEW RADIO OBSERVATIONS 

Short test observations were carried out with the ATCA in the con- 
tinuum at 3 and 6 cm in December 2003 to determine whether emis- 
sion could be detected. The ATCA has six 22-m diameter antennas 
with five of the antennas located on an east-west track of length 
3 km. The sixth antenna lies 3 km further west giving a maxi- 
mum baseline of 6 km. Further continuum observations were made 
at 3, 6, 13 and 20 cm in 2004-2005. Only OH and SiO masers 
were known for OH 354.88-0.54. Therefore, we investigated the 
H2O maser line at 22235 MHz in November 2004 with a short 
observation that yielded a solid detection. We later followed up 
with a 12-h observation in June 2005 with the intention of imag- 
ing the water maser emission and to search for continuum emission 
at 22235 MHz. 

Table 1 summarizes our new continuum observations. We re- 
fer to the observing bands as 1.35, 3, 6, 13 and 20 cm. The actual 
frequencies used were: 22235, 8640, 4800, 2368, and 1384 MHz. 
We used a phase centre of il h 35 m 02.73 s -33° 33' 19.71" 
(J2000) for the continuum observations, and ll h 35 m 02.73 s -33° 
33' 29.41" for the maser observations, corresponding to the stel- 
lar position. To correct for atmospheric amplitude and phase vari- 
ations, the observations of the source were interleaved with short 
observations of a secondary calibrator source (1729-37 for 3, 6, 13 
and 20 cm, and 1710-269 for 1.35 cm). The flux density scale and 
bandpass calibrations were applied from observations of the source 
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Table 1. Log of ATCA observations. 



Date 


Project 


Config 


Time 


Bands 
cm 


Cont/line 


Bandwidth 
MHz 


No. 


08 Jan 2004 


CX053 


6A 


6h 


3,6 


Continuum 


128 


32 


14 Nov 2004 


C1339 


750C 


~45min 


1.35 


H2O maser 


32 


512 


21 Dec 2004 


CI 339 


1.5D 


12h 


3,6 


Continuum 


128 


32 


24 March 2005 


CI 339 


6A 


12h 


13,20 


Continuum 


128 


32 


21 June 2005 


C1339 


6B 


~12h 


1.35 


H2O maser 


16 


512 
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Figure 1. 3-cm ATCA continuum image of OH 354.88-0.54. The large cir- 
cle indicates the approximate optical extent of the PN. The ring is 39" in 
diameter, centred on the stellar position which is marked by the small cross. 
Radio contours plotted correspond to -4, -2, 2, 4, 8, 12, 16, 20 a. lcr is 
0.04 mJy beam , The beam size is shown in the bottom left corner. 



Figure 2. 6-cm ATCA continuum image of OH 354.88-0.54 as for Fig.[T] 
Radio contours plotted correspond to -4, -2, 2, 4, 8, 12, 16, 20, 30, 40, 50, 
60 a. lcr is 0.03 mJy beam . 



1934-638 which is the primary flux calibrator for the ATCA. The 
resolutions achieved are given in column(4) ot Table 2. 

o 
o 

The 3, 6, 13 and 20-cm continuum images were made using | 

u 

natural weighting. All 3-cm data were used because there were no a 
indications of broad scale Galactic emission. We cut out 6-cm data 
with uv-distances <2 kA to eliminate a low-level extended wash 
of Galactic emission and this radically improved the imagery. For 
the 13 and 20-cm data we excluded data with uv-distances <3 kA. 
FiguresQ|2||3|and|4]show the images with the position of the star 
and the approximate location of the outer optical edge of the PN 
also plotted. 

Two obvious sources are detected within the PN at 3, 6 and 
13 cm. The south western source (source A) is the stronger and is 
partially resolved with the resolution of the ATCA. The source to 
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Figure 3. 13-cm ATCA continuum image of OH 354.88-0.54 as for Fig.0 
Radio contours have been chosen to emphasize a possible weak structure 
linking the two continuum components and correspond to -1.5, 1.5, 2, 3, 6, 
8, 12, 14,21,30, 35 and 40 mJy beam" 1 , lcr is 0.14 mJy beam" 1 . 
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Figure 4. 20-cm ATCA continuum image of OH 354.88-0.54 as for Fig.Q 
Radio contours plotted correspond to -2, 2, 4, 8, 12, 16, 20<r. lcr is 
0.18 mJy beam -1 . 



the north of the stellar position (source B) is unresolved and signif- 
icantly weaker, with no detection at 20 cm. 

Table 2 presents the flux densities, positions, sizes and off- 
sets from the star, for source A and source B at each wavelength. 
For source A, the size was estimated for each waveband by fitting 
a Gaussian and then deconvolving to remove the ATCA restoring 
beam. For source B, the Gaussian fitting was consistent in all cases 
with a point source. 

A description of Table 2 follows: col.(l) gives the wavelength 
in cm; col. (2) the source; col. (3) the integrated flux density and 
its uncertainty; col. (4) the semiaxes and p. a. of the major axis for 
the beams used to deconvolve the images; col. (5) the same param- 
eters as in col. (4) for the deconvolved sources; col. (6) the best-fit 
positions to the Gaussians; and col. (7) the modulus of the vector 
distance between each source and the central star. Our upper limits 
are 5a values. 

Source A has a definite nonthermal spectrum with the detected 
flux density falling with increasing frequency. The 20-cm detection 
of 7.7 mJy is comparable to the 10.8±0.6 mJy detection of the 
MGPS1 object (Fig.[3J. Therefore, we associate A with the 843- 
MHz (35-cm) source. In Figure 5, inverse-variance weighted least 
squares fits to power laws are plotted for sources A and B. Both 
sources have strongly nonthermal spectra (F„ oc v mdex ) with in- 



dices of — 0.81±0.01 and — 0.95±0. 11 respectively, probably due 
to synchroton emission. Fitting just the four ATCA data points for 
source A still yields a slope of — 0.81±0.01, indistinguishable from 
the power law fit including the MGPS1 datum, vindicating our 
identification of source A as the 843-MHz source. The power law 
fitted to A is also consistent with the VLA 2-cm upper limit ob- 
tained by Zijlstra et al. (1989) and our own limit on continuum at 
1.35 cm (see below). 

Source B is much weaker than A, but the slope of the radio 
continuum emission is defined by robust detections at 6 and 13 cm. 
Source B appears to suffer substantial absorption at 20 cm (Fig.|5J 
which we interpret as free-free absorption. This absorption is likely 
to occur within the volume of the PN suggesting that source A is 
seen on the front of the nebula while source B is seen at the back. 

The 3-cm flux density of 0.60 mJy detected for source B is 
somewhat higher than that expected for synchroton emission alone 
(0.31 mJy). (We did not include this point in the derivation of the 
spectral index for B because of the high significance of the devi- 
ation from the slope.) We interpret this detection as indicating the 
likely presence of some additional weak thermal radio emission 
(§5). 

In Figure 3, a weak structure appears to provide a link between 
sources A and B at 13 cm (Fig.0. However, this is only a 2a de- 
tection and is not confirmed. There is no evidence for diffuse radio 
continuum emission filling the entire PN. 

The projected mean distances from the stellar position to the 
peaks of sources A and B are the same, within the uncertainties: 
10.4±0.7" (A) and 9.9±0.8" (B). 

Figure|6|represents the discovery of H2O maser emission from 
OH 354.88-0.54, from a 40-minute observation taken in Novem- 
ber 2004. Two peaks are seen at velocities (LSR) of —4.8 and 
+22.2 km s - , significantly displaced by 1-2 km s~ from the cor- 
responding OH maser peaks (—5.5 and +24.3 km s - , respectively: 
CPC). The blue peak is the stronger with a flux density of 0.31 Jy; 
the red peak has a flux density of 0.21 Jy. From these limited data 
we were unable to measure the maser positions. We can infer only 
that we observed two blobs of water maser emission, but we cannot 
interpret the peak velocity differences in terms of relative locations 
of the masing species. It is also noteworthy that the water maser 
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Table 2. Summary of best-fit continuum results to single Gaussians for sources A and B. 



Band 




Tf~\t Q 1 fl 1 1 Y — 1— TTTl G 
lUldl 1111. AZIZ 11 11 3 


IVCSLUllllg UCd.111 


LJC\-.\Jll\ (Jl VCU 


111 


Offset 


cm 




mjy 


semi-axes 


size 


Position 


from star 








(arcsec)(deg) 


(arcsec)(deg) 




arcsec 


1.35 




<1.65 


1.1x0.6 




No continuum detected 




3 


A 


1.86±0.04 


3.0xl.83(-20.9) 


4.48xl.39(-12.8) 


17:35:01.74 -33:33:36.0 


10.5±1.4 


6 


A 


2.98±0.03 


5.1x2.7(-15.7) 


4.4x0.84(-12.5) 


17:35:01.74 -33:33:35.9 


10.5±1.4 


13 


A 


6.02±0.14 


9.3x4.5(2.3) 


6.0xl.9(-1.6) 


17:35:01.75 -33:33:35.6 


10.2±1.4 


20 


A 


7.68±0.18 


13.9x7.3(0.9) 


2.8x0.65(-41) 


17:35:01.74 -33:33:36.0 


10.5±1.4 


3 


B 


0.60±0.04 


3.0xl.83(-20.9) 


Point source 


17:35:02.73 -33:33:19.3 


10.1±1.4 


6 


B 


0.54±0.03 


5.1x2.7(-15.7) 


Point source 


17:35:02.75 -33:33:19.5 


10.0±1.4 


13 


B 


1.06±0.14 


9.3x4.5(2.3) 


Point source 


17:35:02.64 -33:33:19.7 


9.7±L4 


20 


B 


<0.90 


13.9x7.3(0.9) 


Point source 


No continuum detected 
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Figure 5. Spectral energy distributions of sources A and B fitted by power 
laws. Filled triangles indicate upper limits at 2 cm from the VLA and 1.3 cm 
from our data for A, and at 20 cm, from our data, for source B. 



profile is significantly above zero in the centre, similar to the OH 
profile (CPC). 

Despite the 12-h track obtained during the June 2005 ATCA 
observation, we had a non-detection of any water maser emission 
with a 5<7 limit of ~40 mJy. Extrapolating from Figure 3 of CPC 
we estimate that we detected water masers at an OH maser phase 
of 0.57 but failed at a phase of 0.72. This interval is such a small 
fraction of a period that the source must be highly variable and 
the water masers clearly do not follow the very regular OH maser 
variations. The range of variation in H2O maser intensity from the 
2004 to the 2005 observation was about a factor of 8. We looked 
for radio continuum emission at 1.35 cm using the spectral line data 
but achieved no detection, with a 5<r upper limit of 1.65 mjy. 




Figure 6. Spectrum of H2O maser emission detected from the PN. 
4 OPTICAL IMAGERY 

Figure represents a quotient map obtained directly from Super- 
COSMOS data for the PN by dividing the Ha pixel data by the cor- 
responding data from the short-red (SR) exposure. However, it also 
involves use of the Bond et al. (2002) technique of variable point- 
spread-function (PSF) matching between the two images prior to 
division. This enhances the faint nebula by effectively removing 
most of the emission from point sources in the field. The figure 
shows that the outer boundary of the PN is quite sharply defined 
around most of its circumference. Emission can clearly be seen 
above the sky over most of its projected area although the bright- 
ening of the south-eastern region is particularly striking. There is 
also the suggestion of a northward extension although this is very 
faint and merges imperceptibly with the interstellar medium (ISM). 
The north-south extent could be as much as 20% greater than the 
east-west diameter in this newly processed image. Fig.|8|shows the 
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Right Ascension (J2000) 

Figure 7. Ha/SR quotient image of the PN with variable PSF-matching 
before dividing the Ha and SR data. The linear feature east of the PN is the 
track of a satellite. 



relationship between sources A and B and the PN in the plane of 
the sky. 

An old UKST Ha plate which includes the PN region exists 
in the archive. This exposure (HA2307) was taken in 1976 for John 
Meaburn on fast but coarse-grained O-9804 emulsion with a mo- 
saic Ha filter. The enhanced south-east section of the PN is clearly 
visible on this old exposure but the remaining boundary is hard 
to discern. This, coupled with the poorer resolution of this expo- 
sure, make it essentialy impossible to measure whether there has 
been any expansion of the nebula ring over the intervening 20-plus 
years. 



5 DISCUSSION 

Nonthermal emission from early-type stars was found in Wolf- 
Rayet stars (WRs) by Abbott, Bieging, & Churchwell 1984; Ab- 
bott et al. 1986; Bieging et al. 1989). Its origin lies in the acceler- 
ation, by strong shocks, of a small population of relativistic elec- 
trons (White 1985). In WR binary systems these shocks arise from 
the collision of the separate winds of the two components. It is 
less certain whether the phenomenon occurs in single WR stars. 
The paradigm applied to OH 354.88-0.54 involves a single star but 
two stellar winds from different phases of evolution. Absorption of 
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Figure 8. An enlarged version of Fig.lTloverlaid by white contours of 6-cm 
radio continuum emission corresponding to 0.3, 0.5, 0.8, 1.1, 1.4, 1.6, 1.9, 
2.2 mjy beam -1 . White circle marks the approximate outer boundary of 
the PN; cross marks the position of the star and OH masers. 



synchrotron emission has been observed in colliding winds from 
WRs by Chapman et al. (1999), who attributed this to thermal ab- 
sorption by gas along the line-of-sight through the ionized wind. 
The observed continuum energy distribution for WR 48 (Chapman 
et al. 1999, Fig. 11) reveals nonthermal emission between 3 and 
13 cm but only an upper limit at 20 cm, implying a pronounced 
turnover longward of ~13 cm. These authors separated the ther- 
mal and nonthermal contributions to the spectrum of WR48 and 
found the nonthermal component to have a slope of —0.83. This 
slope matches that for OH 354.88-0.54 sources A and B (within 
the errors), strongly suggesting that we may be viewing nonthermal 
emission in this PN that is also caused by the collision of two winds, 
the low density, high velocity PN fast wind and the high density, 
slowly moving original PN ejecta. Synchrotron self-absorption is 
unimportant in sources with brightness temperature below 10 10 K 
(Condon 1992) but thermal absorption within the PN is still relevant 
and could cause the abrupt turn-down in the spectrum of source B. 
Locating source A near the front side of the PN, as we view it, 
would account for the lack of absorption in its radio spectrum. 



5.1 Estimates of electron density 

The PSF-subtracted image of the PN (Fig.0 indicates that the Ha 
brightness varies significantly across the PN suggesting that the 
electron density also varies. It is also possible that the southern rim 
of the PN is brightest because the PN's slow wind suffers a shock at 
the interface with the ISM. Source A lies within the western portion 
of the brightest (southern) region of Ha emission, for which CPC 
determined N e ~4000 cm -3 , whereas B appears in an area of much 
lower Ha brightness to the north (Fig.|7J. We can use the turnover 
in source B's radio spectrum to estimate the mean electron den- 
sity through the nebula in its direction. The thermal optical depth 
through the nebula may be calculated following Mezger & Hen- 
derson (1967, eqn (5)), and rearranged into an expression for the 
turnover frequency (in GHz), ^t~0.3(T- 1:35 N^L) ' 5 (Lang 1999, 
eqn. 1.224), where L is the path length through the PN in pc (we 
adopt the diameter of the PN, 1.9xl0 18 cm, 0.6 pc, for L: note 
that CPC's Table 1 gave the diameter as 0.3 pc; this was in fact 
the radius). The turnover wavelength must lie beyond 13 cm but be 
less than 20 cm, so vt must lie between 1.5 and 2.3 GHz. Taking 
T e ~10 4 K, we conclude that the mean nebular N e along the line-of- 
sight toward B is between 3200 and 4900 cm -3 , in excellent agree- 
ment with the density derived from the [Sll] lines (4000 cm -3 : 
CPC). 



5.2 Possible free-free emission in source B 

We also consider whether the apparent excess of 3-cm emission in 
B, above the nonthermal slope, might be due to free-free emission, 
as seen in WR stars. Firstly we have separately analyzed the two 3- 
cm data sets to investigate the reality of this excess by fitting a point 
source to the calibrated u-v data. This should give a more reliable 
flux estimate than fitting in the image plane, if there were no con- 
fusion and the emission were unresolved. We find 0.61 ±0.03 mjy 
(Jan 2004) and 0.59±0.03 mJy (Dec 2004), where the errors are 
the rms residuals from the fitting process alone. Our measurement 
of the same value at two epochs separated by nearly a year gives 
strong support for the reality of the measured flux density. Extrap- 
olating source B's nonthermal spectrum to 3 cm implies a level of 
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0.31 mjy. Therefore, we conclude that B shows a real excess at 
3 cm of 0.29 mjy. 

To estimate the expected radio free-free level we have ap- 
plied an absolute calibration to the Ha image of the nebula ac- 
cording the method described by Pierce (2006), in which a scaled 
continuum (SR) image is subtracted from the Ha image, and a 
factor is applied to relate the SHS difference image to the spa- 
tially much lower-resolution, but absolutely calibrated, Southern 
Ha Sky Survey Atlas (SHASSA: Gaustad et al. 2001). For this 
PN, the SR scale is 0.99 and the conversion factor in this field 
is 14.9 counts pixel -1 per Rayleigh. Integrating the entire PN 
spatially, subtracting background sky (assessed as the average of 
that measured in four different locations to provide a robust back- 
ground) and three point sources interior to, or on the rim of, the 
nebula, and after correction for the 36% contribution of the red 
[Nil] line to the emission in the bandpass of the SHS Ha fil- 
ter, we derive 3.6 Rayleighs in Ha alone. This is equivalent to 
2.4xl0~ 14 erg cm -2 s~ 1 . We now require the total extinction 
to the PN. Based on Fitzgerald (1968), Lucke (1978) and espe- 
cially on Neckel & Klare (1980), we find Ay«3.5 in this direc- 
tion, corresponding to Aj/a ~2.8. Thus the intrinsic Ha flux is 
3.3 xl0~ 13 erg cm -2 s -1 , and the corresponding H/3 flux would 
be l.lxlO" 13 erg cm" 2 s" 1 . 

This is convertible into the radio free-free flux density using 
Condon (1992: eqns. (3) and (4a)). The estimated free-free emis- 
sion at 3 cm, for the entire nebula, is 0.33 mjy and this is compa- 
rable to the excess emission detected at 3 cm. Source B was un- 
resolved and thus only occupies a small part (less than about 2%) 
of the nebular area. For a uniform distribution the predicted free- 
free emission from source B would be less than 0.01 mjy. How- 
ever, one would expect the radio thermal emission to follow the Ha 
brightness distribution and, therefore, to be patchy. There is prece- 
dent for this in the Luminous Blue Variable HD 168625 (Leitherer 
et al. 1995), where the radio continuum shell closely mimics the 
morphology of the partially filled Ha shell. If the free-free emis- 
sion were highly clumped toward source B, then a larger fraction 
of the total thermal emission might be concentrated there. Alter- 
natively, if our line-of-sight to source B through the PN suffered 
significant internal extinction then the estimated intrinsic Ha emis- 
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sion would also increase. Consequently, local overdensities of ion- 
ized gas could be offset by internal nebular extinction due to dusty 
clumps. For example, if 20% of the total observed free-free emis- 
sion were located within source B (an overdensity greater than a 
factor of 10), and there were an extra extinction (Ay) of 2.1 mag- 
nitudes (Ah a of 1.7 mag) within the nebula, then radio free-free 
emission of 0.3 mjy would still be observed. (In this case the to- 
tal PN free-free emission would be 1.65 mjy, boosting the intrinsic 
Balmer line fluxes by x5, requiring a line-of-sight Ahc of 4.6, 
equivalent to an Ay of 5.6 mag, of which 3.5 mag arise from the 
ISM and 2. 1 mag within the PN. 

The upper limit on B at 20 cm also excludes a cut-off in the 
electron energy spectrum because there is no intermediate regime 
with S^oci/ 1 / 3 between the optically thin and thick spectral indices 
(Ginzburg & Syrovatskii 1965). The fact that A and B have almost 
the same nonthermal slope strengthens the case for either a com- 
mon volume in which we observe this nonthermal emission, or two 
separated volumes in which very similar strong shocks generate 
relativistic electrons. The fast wind, therefore, must encounter es- 
sentially the same physical conditions at the two locations in which 
it impacts the nebular ejecta. 



5.3 Background source contamination 

Both sources A and B are located within the ring and no other 
strong background objects are apparent. One can ask whether 
sources A and B are truly associated with the PN or might represent 
a double radio source viewed by chance along the line-of-sight. In- 
spection of our observed images covering the ATCA primary beams 
at 20 cm (~1000 arcmin 2 ) and 6 cm (~80 arcmin 2 ) suggests that 
7 and 3 sources, respectively, are recognizable above the noise, at 
all flux densities, and none appear double. Calculating the formal 
extragalactic confusion from Log N— Log S relationships at 20 and 
6 cm (Ledden et al. 1980; Bridle 1989) one would expect to detect 
15 radio sources within the ATCA primary beam at 20 cm with a 
flux density greater than the sum of A and B, but less than one simi- 
larly within the 6-cm primary beam. The probabilities of even a sin- 
gle background source as bright as A appearing within the PN Ha 
image are approximately 1/500 (6 cm) and 1/200 (20 cm). There 




=J I I I I I I I I I I I I I I I - 

1 7 n 35 m 0Z.G a ai.s a Ol.fl 3 Dlji a 

Right Ascsnnion (JEOOO) 

Figure 9. 6-cm continuum emission of source A shown both as a greyscale 
image and overlaid by its own contours in white, for levels of 0.2, 0.4, 0.7, 
0.9, 1.2, 1.4, 1.6, 1.85, and 2.0 mjy beam- 1 . The black dashes mark the 
locus of the ends of the major axes of the white contours. Note the distinct 
curvature of this locus. 



are several other arguments that imply a physical linkage with the 
PN: they have statistically the same nonthermal radio spectral in- 
dices instead of the more common thermal sources found along the 
Galactic plane; there may be a faint connection between the emis- 
sion of A and B; and A lies very close to the edge of the PN and is 
clearly curved, roughly concentrically, with the rim (Fig.|9j, mak- 
ing a physical interaction between A and the ejecta that define the 
PN rim highly plausible. 

5.4 Radio luminosity and kinetic energy: comparison with 
WR stars 

We have compared our continuum observations at different epochs 
and see no indications of variability at any of the wavelengths ob- 
served over a period of 1.5 yr. Chapman et al. (1999: their Table 3) 
tabulated the monochromatic 6-cm luminosity of WRs that emit 
nonthermal radio emission. The average of the 9 WRs detected 
was 2.1±0.6xl0 19 ergs s" 1 Hz" 1 . The same calculation yields 
3.7 and 0.7 x 10 19 ergs s -1 Hz -1 for A and B, respectively. Sum- 
ming the two nonthermal sources in OH 354.88-0.54, the PN is 
more luminous than the average nonthermal WR emitter. Adding 
A and B we find that the approximate nonthermal radio luminosity 



(for an adopted bandwidth of 10 GHz; cf. Chapman et al. 1999) is 
4.4xl0 29 ergs s" 1 Hz" 1 . High energy emission is detected from 
some WRs in which colliding winds occur. OH 354.88-0.54 has not 
been detected based on a search of X-ray and 7-ray catalogues. In 
particular, Chandra has not observed in its direction. 

We now compare OH 354.88-0.54 with WR stars, for which 
a typical radio luminosity is a few xlO 29 ergs s" 1 . The total 
kinetic energy (KE) flux calculated from the known mass loss 
rates (M) of WRs (a few x 10~ 5 M yr _1 ) and wind velocities 
(V win d~2000 km s" 1 ), \ M W wind 2 , is typically 6x 10 37 ergs s" 1 . 
For a WR— O-star collision about 0.3% of the flux passes through 
the interaction region, so the KE flux at the interaction is about 
2x 10 35 ergs s" 1 (see Chapman et al. (1999) for details). 

For the detected continuum sources A and B, the 5-GHz ra- 
dio luminosity is 4.4x 10 19 ergs s _1 Hz -1 and the total integrated 
radio luminosity (h ra dio) is thus about 4.4x 10 29 ergs s -1 . We as- 
sume that the ratio of KE flux/radio luminosity is the same as for 
WR stars (about 10 6 :1). Thus approximately 4.4 xlO 35 ergs s _1 
should be available at the shock regions. By comparison with WR 
binaries, the total KE flux of the fast wind of OH 354.88-0.54 is not 
well determined as we do not know what fraction of the fast wind 
encounters the shocked regions. Source A occupies about 1/100 of 
the total shell area (estimated at 6 cm) while source B is smaller. 
For a smooth distribution this suggests that the KE could be about 
lOOx higher than at the shock. However, if the fast wind is clumpy 
or asymmetric then this is likely to be considerably overestimated. 

As an example, we consider a scenario where 10% of the total 
KE flux reaches the detected shocked regions giving a total KE flux 
of 4.4x 10 36 ergs s _1 . For V/ aa t, the velocity of the fast wind, we 
adopt 1000 km s _1 . From these values we can crudely estimate the 
mass-loss rate of the fast wind to be M/ as t«1.2x 10 M0 yr _1 . 

We now estimate the velocity of the compressed shell (V shell) 
for this mass-loss rate. From eqn.(9) of Zijlstra et al. (2001) and 
using M agb = 7xl0~ 5 , W agb = 15 km s" 1 , M /ost = 1.2xl0- 5 , 
and V fast = 1000 km s , we have /i = 0.17 and £ = 67, where 
ji = Mf ast / M-agb and £ = V/ as t/V asi) . For these values the velocity 
of the compressed shell is ~60 km s~ and this can be compared 
to the observed line-of-sight optical velocity of 20 km s _1 . If less 
than 10% of the total energy flux reaches the shocked regions then 
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both Mf a st and V shell would be higher while, if a larger fraction 
of the energy flux were intercepted, then Mf ast and V s h e u would 
be lower. 

5.5 Wind interactions 

Despite the roughness of these estimates they do provide an idea of 
the relevant interactions. The fast wind ploughs into the slow AGB 
wind creating shocks, non-thermal emission and a compressed shell 
on the inner face. If the fast wind has a momentum (M f ast X V/ os t) 
that is roughly comparable to (or greater than) that of the slow wind 
then, when they collide, the slow wind is compressed and acceler- 
ated and this increase in velocity is indicated by the velocity of 
the optical emission. The shock excitation must be stronger in the 
compressed intranebular shell than at the ISM and the rather high 
density derived from the red [Sll] doublet is more suggestive of a 
post-shocked region rather than a typical PN envelope. Therefore, 
it is more likely that the optical emission in the PN arises at the in- 
terface of the fast wind and AGB wind, rather than at the interface 
of the AGB wind and the ISM. 

The geometry of the emerging fast wind is unknown because 
any spherically symmetric flow from the star would be immedi- 
ately modified by its encounter with the dense circumstellar OH 
toroid. Naively, were this toroid to constrain the fast wind to the po- 
lar directions (north-south) in which OH maser emission is weak- 
est (Welty, Fix & Mutel 1987), one might expect to see nonther- 
mal "hot spots" in two opposing directions. There is evidence of a 
10-fim source around the star that is significantly elongated north- 
south (Cobb & Fix 1987) but this need not trace any dynamical 
structure. In fact, Figure Q suggests that sources A and B do not 
correspond to a bipolar wind encountering a spherical distribution 
of slowly-moving ejecta. Consequently, one might have to invoke 
either a fast wind capable of breaking out of the OH toroid in any 
region of locally reduced density, or a non-spherical distribution 
of the original PN ejecta, or both. However, Table 2 shows that 
both sources appear projected in the sky at the same radial distance 
from the central star. This is consistent with a spherical ejection 
for the slow wind. The natural direction for any asymmetry of the 
fast wind or of the original debris would be in the direction roughly 
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from northwest to southeast of the star, based on the morphology 
of the Ha nebular brightness and the line bisecting sources A and 
B through the possible weak bridge of emission observed at 13 cm 
(Fig-EJ. There is also evidence for a preferred axis in the same p. a. 
in the OH velocity slices between —5 and —8 km s _1 published 
by Welty et al. (1987: their Fig. 15). At more negative velocities 
the OH brightness indicates weaker OH emission along a north- 
south axis. None of these directions appears to be associated with 
the nonthermal shock zones. 

Strong nonthermal radio continuum emission has also been 
observed from the OH/IR star OH 326.53-00.419 (also known as 
D046), with a spectral index of —0.8. The OH maser emission from 
D046 occurs over an unusually broad range of 80 km s _1 (Seven- 
ster & Chapman 1999; Sevenster & Chapman 2005). D046 appears 
to represent a post-AGB object, in transition between the AGB and 
PN phases. These nonthermal emitters may all arise because of the 
impact of fast winds on circumstellar debris near the end of AGB 
life. OH 354.88-0.54's intermediate mass progenitor may confer its 
uniqueness, in that it is still a pulsating long-period variable after it 
has already created a visible PN, simply because the time scale for 
stellar evolution is much shorter for a 4 Mq star than for a solar- 
mass object. 

6 CONCLUSIONS 

The ionization of the PN associated with OH 354.88-0.54 appears 
to arise through shock interactions between the AGB star's fast 
wind and previously ejected cold nebular material. Our discovery 
of patches of nonthermal radio emission in this PN suggests a kin- 
ship with colliding-wind WR binaries but, in OH 354.88-0.54, the 
two winds were produced by a single central star. The locations of 
the shocked, nonthermally emitting, regions do not suggest a spher- 
ical interaction zone between the winds, nor do they align with any 
geometrically meaningful direction in the PN itself. However, the 
nonthermal regions each show the same projected distance from 
the central star in the plane of the sky, and there is the possibility 
that faint radio emission links the two regions. These characteristics 
suggest that the anisotropic fast wind collides with the continuous 
inner spherical surface of the slow wind. Therefore, this PN offers 



a rare opportunity to recognize the nonspherical geometry of the 
fast wind, close to the time of its onset. 
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